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2\ 



We show analytically that the QCD potential can be expressed, up to an 0(Aq CD 
uncertainty, as the sum of a "Coulomb" potential (with log corrections at short distances) 
and a linear potential, within an approximation based on perturbative expansion in as 
and the renormalon dominance picture. The expansion of Vqcd(?") is truncated at 0(a>g ) 
[iV = 67r/(f3oas)], where the term becomes minimal according to the estimate by NLO 
renormalon, and is studied for N S> 1. Analytic expressions for the linear potential are 
obtained in some cases. 



1 Introduction 



Analyses of the static QCD potential Vqcd(^) within perturbative QCD entered a new phase 
when the cancellation of the leading order (LO) renormalons between the QCD potential and 
the pole masses of quark and antiquark was discovered [1]. Convergence of the perturbative 
series improved dramatically and much more accurate perturbative predictions became avail- 
able. Subsequently, several studies [2-6] showed that perturbative predictions for Vqcd(^) agree 
well with phonomenological potentials (determined from heavy quarkonium spectroscopy) and 
lattice calculations of Vqcd(?")> once the LO renormalon contained in the QCD potential is 
cancelled (see also [7]). In fact the agreement holds within the perturbative uncertainty of 
(9(AQ CD r 2 ) estimated from the residual next-to-leading order (NLO) renormalon [8]. Despite 
of different prescriptions used for cancelling the LO renormalon, all these perturbative predic- 
tions were mutually consistent within the 0(Aq CD r 2 ) uncertainty* These observations indicate 
validity of the renormalon dominance picture for the QCD potential. 

Empirically it is known that phenomeno logical potentials and lattice computations of Vqcd(?") 
are both approximated well by the sum of a Coulomb potential and a linear potential in the 
range r > 0.5 GeV -1 [9]. The linear behavior at large distances is consistent with the quark 
confinement picture. For this reason, before the discovery of the renormalon cancellation, it 
was often said that perturbative QCD is unable to explain the "Coulomb-plus-linear" behavior 
of the QCD potential. 

Once the cancellation of the LO renormalons is incorporated, the perturbative QCD poten- 
tial gets steeper than the Coulomb potential at large distances. This feature can be understood, 
within perturbative QCD, as an effect of the running of the strong coupling constant [10,2,3]. 
On the other hand, it is not obvious whether the QCD potential is rendered to a "Coulomb- 
plus-linear" form by this effect. The perturbative uncertainty due to the residual renormalon 
is of 0(r 2 ), hence there is a possibility that the 0(r) term of the potential at r < Aqq D is pre- 
dictable within perturbative QCD. In this paper, by considering a certain limit of a finite-order 
perturbative expansion of Vqcd(^) based on the renormalon dominance picture, we show that 
indeed the potential can be decomposed into a "Coulomb-plus-linear" form, up to an C(AQ CD r 2 ) 
uncertainty. Our prescription gives a prediction consistent with the previous predictions [2-6]. 

In Sec. 2 we set up our conventions for our analysis. Sec. 3 presents an analysis in the 
large-/?o approximation; Sec. 4 presents an analysis based on renormalization-group (RG), 
incorporating 1-, 2-, and 3-loop running of the coupling constant. Discussion and conclusions 
are given in Sees. 5 and 6, respectively. 

2 Perturbative QCD potential and renormalons 

The static QCD potential is defined from an expectation value of the Wilson loop as 



*This is true only in the range of r where the respective perturbative predictions are stable, since all the 
perturbative predictions go out of control beyond certain distances. 
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where V is a rectangular loop of spatial extent r and time extent T. The second line defines 
the l^-scheme coupling contant, av(q), in momentum space; Cf = 4/3 is the second Casimir 
operator of the fundamental representation. In perturbative QCD, ay(q) is calculable in a 
series expansion of the strong coupling constant: 
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a v (q) = a s £ P n Qn(jji/q)) (^-J = a s (q) f] a n (^f^)" 5 «n = ^(0). (2) 

Here, P„(£) denotes an n-th-degree polynomial of £. In this paper, unless the argument is 
specified explicitly, as = as(/f) denotes the strong coupling constant renormalized at the 
renormalization scale /i, defined in the MS scheme. The series expansion of as{q) in terms of 
as is determined by the RG equation 

q ^ a 5 (,) = E^(-^J , (3) 

where /3 n represents the (ra + l)-loop coefficient of the beta function. Therefore, at each order of 
the expansion of ay(q) in as, the only part of the polynomial P n (ln(/x/ g)) that is not determined 
by the RG equation is a n . The above equations fix our conventions. 

It is known [11] that a n for n > 3 contain infrared (IR) divergences. We will discuss this 
issue in Sec. 5, whereas in Sees. 3 and 4 we treat a n as finite constants. 

According to the renormalon dominance picture, the leading behavior of the 0(a^ +1 ) term 
of VqcdM at large orders is given by the LO renormalon contribution as Vqc D (r) ~ const, x 
n\ ((3 as/(2ir)) n n 5 / 2 , where S — Pi/Po [12]. In the computation of the heavy quarkonium spec- 
trum, the LO renormalon gets cancelled against the LO renormalons contained in the quark and 
antiquark pole masses. Considering this application, if we subtract the LO renormalon contribu- 
tion from VqQ B (r), its large-order behavior becomes const, xr 2 n ! ( / 3 a5/(67r)) ri n 3 ' 5 / 2 due to the 

NLO renormalon contribution. Then Vqq D (?') (after the LO renormalon is subtracted) becomes 
minimal at n ~ iV = 6rr/ (Po a s) an d its size scarcely changes for iV — <C n <^ N + \^N. 

In view of the usual property of asymptotic series, we simply truncate the series expansion 
of the potential at the order where the term becomes minimal according to the renormalon 
dominance picture, i.e. at 0(ag): 

V N (r) = [V Q cb(t)] n = -47TCV / (0 ^ M?)]n- (4) 

Here and hereafter, [X]n denotes the series expansion of X in 0:5 truncated at 0(a^). The 
purpose of this paper is to examine VJv(r) for ^> 1 while keeping [13] finite, using certain 
estimates for the all order terms in Eq. (2). The motivation for considering the large N limit is 
that it corresponds to the limits where the perturbative expansion becomes well-behaved (small 
expansion parameter) and where the estimate of Vq^ D (r) by renormalon contribution becomes 
a better approximation around n ~ N. Note that large N corresponds to small as and large \x 
due to the relation between N and as. 
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Clearly, Vn{t) cannot be written as a "Coulomb-plus-linear" form for finite N, since it is 
given as the Coulomb potential (— Cp<ys/ r ) times an (iV — l)-th-degree polynomial of ln(/xr), 
and therefore, Vjy(r) — > as r — > oo. We will see, however, that Vjv(r) tends to a "Coulomb- 
plus-linear" potential (plus a quadratic potential) in the large N limit. 



3 Vjy(r) in large— j3 approximation 

The large-/9 approximation [14] is an empirically successful method for estimating higher-order 
corrections in perturbative QCD calculations. For Vqcd(?")j this approximation corresponds to 
setting a n = (5/3o/3) n and all /3 n = except /3q. (Therefore, it includes only the one-loop 
running of as(q)-) In this section, with these estimates of the all-order terms of Vqcd(^), 
we examine Vjv(r) for N ^> 1. The reasons for examining the large-/3o approximation are as 
follows. First, because this approximation leads to the renormalon dominance picture; in fact, 
the renormalon dominance picture has often been discussed in this approximation. Secondly, 
as stated in Sec. 1, the running of the strong coupling constant makes the potential steeper 
at large distances as compared to the Coulomb potential; hence, we would like to see if the 
potential can be written as a "Coulomb-plus-linear" potential when only the one-loop running 
is incorporated as a simplest case. We first present the results, discuss some properties, and 
then sketch how we derived our results. 

Results 

We define A = e 5/6 A^| op , where A^| op = nexp[-2n/(p a s )]. In this section, we assume 
e~ 5/6 iT 1 = A -1 e -JV / 3 < r < A~ 1 e N/3 when taking various limits. Note that, as A" — > oo, the 
lower bound (e _5 / 6 /i _1 ) and the upper bound (A~ 1 e N ^ 3 ) of r go to and oo, respectively. 

Vn(t) for N ^> 1 within the large-/3o approximation can be decomposed into four parts 
corresponding to {r _1 , r°, r 1 , r 2 } terms (with logarithmic corrections in the r^ 1 and r 2 terms): 

V<*\r) = ^Av(Ar,N), (5) 
v(p, N) = v c (p) + B(N) + Cp + D(p, N) + (terms that vanish as N -> oo). (6) 
(i) "Coulomb" part: 



7T 1 f°° 

vc{p) = 1 — / dxe x arctan 

p p Jo 



(7) 



" vr/2 
_\n(p/x) 

where arctan a; G [0, 7r). The asymptotic forms are given by 

( Vci f[ ~ " " -* ° (8) 

[ «c(p) ~ - p, p ^ oo 

and both asymptotic forms are smoothly interpolated in the intermediate region, 
(ii) constant part''": 



f°° r / 3 \ ./v 



^The 0(1/N) and 0(1/N 2 ) terms in eq. (9) are irrelevant for N — > oo. We keep these terms in B(N) for 
convenience in examining vjf°\r) at finite iV; see Fig. 1 below. 
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(iii) linear part: 
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(iv) quadratic part: 
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where 9(x) is the unit step function and ^ E = 0.5772... is the Euler constant. The asymptotic 
forms of d(p) are given by 



d (p) ~ [21nlni + ln| + 7E 
~ [21nlnp + ln| + 7 jB 



P 
P 
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(13) 



and in the intermediate region both asymptotic forms are smoothly interpolated. 

Although the constant part of Vjf°\r) diverges rapidly as N — > oo, the divergence can 
be absorbed into the quark masses in the computation of the heavy quarkonium spectrum. 
Therefore, in our analysis, we will not be concerned with the constant part of the potential but 
only with the r-dependent terms. 

The quadratic part of vjf°\r) diverges slowly as A 3 r 2 lniV ~ A 3 r 2 \n\n(p/A)} We may 
consider this feature to be a characteristic property of renormalons for the following reasons. 
(1) If the series expansion of m po \ e (m^, as) or Vqcd(^) is truncated at the order corresponding 
to the minimal term of^the LO renormalon contribution, i.e. N' = 2n /(j3 as), [mpoiejiv or 
[Vqcd(^)]w diverges as AlniV' (within the large-/3o approximation). We may compare AlnA^' 
with the usual interpretation that m po i c and Vqcd(^) contain 0(A) perturbative uncertainties 
due to the LO renormalons. (2) We have checked that even if we incorporate the effect of 
the two- loop running, i.e. even if we set j3i ^ 0, the quadratic part of Vjv(r) still diverges as 
A 3 r 2 lnAf. Therefore, we interpret that the quadratic part of vjf°\r) represents an 0(A 3 r 2 ) 
uncertainty, following the standard interpretation on the perturbative uncertainty induced by 
the NLO renormalon. In this respect, we note that the dependence of vjf°\r) on is mild; 
for instance, as shown in Fig. 1, the variation of vjf°\r) is small (after the constant part is 
subtracted) in the range r < A^ 1 as we vary A^ from 10 to 100; it corresponds to a variation of 
p/A^ op from 30 to 3 x 10 14 . 

The "Coulomb" part and the linear part are finite as A^ — > oo. In Fig. 1, we see that 
vjf°\r) is approximated fairly well by the sum of the "Coulomb" part and the linear part (up 

■^Within the potential-NRQCD framework, this divergence or scale-dependence can be absorbed into the 
0(r 2 ) term of a non-local gluon condensate in the operator product expansion [15]. 
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Figure 1: Truncated potential after the constant term is subtracted, v(p,N) — B(N), (dashed black) 
vs. p for N = 10, 30, 100. "Coulomb-plus-linear" potential, vc(p) + Cp, (solid red) is also plotted, 
which is hardly distinguishable from the A^ = 30 curve. 



to an r-independent constant) in the region r < A -1 when we vary N between 10 and 100. 
Moreover, as long as y^lnA^ < 0(1), the difference between vjf°\r) and the "Coulomb-plus- 
linear" potential remains at or below 0(A 3 r 2 ) in the entire range of r. 

Outline of derivation 

Let us write L = ^ a In ( ) = 1 + ^ln^j. vjf°\r) is denned as the Fourier transform 

of (^C F a s /q 2 )En=o Ln = (4nC F a s /q 2 )(l - L N )/(l - L). After integration over angular 
variables, it is given by 



V, 



(A), 



N 



2C F a s f°° , sin(gr) 1 - L N 

/ dq — 

7r jo qr 1 — L 



0o 



Vx(Ar) + f 2 (Ar, N) 



(14) 



where we separated the integral into two parts after deforming the integral contour slightly: 
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(15) 



1+ A7 ln V. 



kp ln(l/fc) 



1 + 



v 



(16) 



i>i and t>2 are defined by the first equalities of (15) and (16), respectively. Contributions from 
the pole at k — 1 in v\ and t>2 cancel, since the original integral (14) does not contain a pole 
at q — A. In the second equality of (15), we deformed the integral contour into the upper half 
plane on the complex fc-plane and integrated by parts. As for t> 2 , since [1 + -j| lnQ)]-^ — ► 1/k 3 
as N — > oo, the constant (p°) and quadratic (p 2 ) terms in the integral become IR divergent in 
this limit. On the other hand, the negative power of k induces the positive power behavior of 
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p, i.e. the linear and quadratic terms, in v 2 in the large N limit. We define 



v 2 (p, N) 



A 



+ B(N) + Cp + D(p, N) + (terms that vanish as N -> oo) 



;i7) 



where -D(p, N) = 0(p 2 ). In computing A and C, it is convenient to first remove divergences by 
subtracting appropriate constant and quadratic terms from v 2 . Let 



v 2 (p, N, k ) 
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ln(l/fc) 
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(18) 



where k = 2 is an IR cutoff. Now we may send iV — > oo before integration over fc. ^(p, oo, fc ) 
is finite for < p < oo and differs from v 2 (p, N) only by constant (p°) and quadratic (p 2 ) terms, 
apart from terms that vanish as N — > oo. One can show that the p _1 and p 1 terms stem only 
from the first term of (18): 



C 



lim pv 2 (p,oo,k ) = -7r, 

, 1 9 2 
hm -7-^7 
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[pv 2 (p, oo,k )] = 
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(19) 
(20) 



The constant and quadratic terms can be calculated directly from v 2 : 

B(N) = \imJ-[pv 2 (p,N)} 
p^o op 



f°° , e~* / 3 \N r-t , e 2t / 9 9 \ 
-/ di — ( 1 + — i ) -/ rft—(l-—t 2 + •••), 

A i V iV / 7-oo t v 2A^ /' 



J D(p,iV)=t; 2 (p,iV) 



- + 5(iV) + Cp 

p 



(21) 
(22) 



In the second equality of (21) we set t = ln(l//c) and expanded the integrand in 1/N in the 
region t < 0. It is then straightforward to obtain (9). One may separate a divergent part as 
N — > oo from (22) in a similar manner. As for the finite part (iV-independent part), we deform 
the integral contour into the upper half /c-plane to obtain (11), (12). Finally i>c(p) is given by 
the sum of A/p and vi(p)} The asymptotic forms of i>c(p) and rf(p) are obtained by expanding 
the integrands in lnx. 

We made a cross check of our results by comparing v(p, N) and Vc{p) + B(N) + Cp+D(p, N) 
numerically for 3 < iV < 1000, after subtracting the divergent terms from both. The difference 
diminishes swiftly with N. 



4 V)y(r) with 1-, 2-, and 3-loop running of cts(q) 

In this section we examine V/v(r) in three cases corresponding to the following estimates of the 
all order terms of VqcdM: 

§Since the leading behavior of vjf°\r) as r — > is const./ (r In r) as determined by the 1-loop RG equation, 
the A/p term of v-i must be cancelled by the 1/p term contained in v\. 
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Figure 2: Integral contours C\ and C 2 on the complex (/-plane, q* denotes the Landau singularity of 
as(q)- For 1-loop running, q* is a pole; for 2- and 3-loop running, q* is a branch point. In the latter 
case, branch cut is on the real axis starting from q* to —00. 



(a) [1-loop running] /3q, a®, exact values, f3 n = a n = (n > 1); 

(b) [2-loop running] (3 , ao, a\\ exact values, (3 n = a n = (n > 2); 

(c) [3-loop running] f3 , /3 2 , a , ai, a 2 : exact values [16,17], /3 n = a n = (n > 3). 

We assume f3 , P±, /3 2 ,a , ai,a 2 (exact) > 0.* In the standard 1-, 2-, and 3-loop RG improve- 
ments of Vqcd(^)) the same all-order terms as in the above cases are resummed; the difference 
of our treatment is that the perturbative expansions are truncated at 0(a<[). We note that the 
renormalon dominance picture is consistent with the above estimates of higher-order terms, or 
more generally, with the RG analysis [12]. All the results for the case (a) can be obtained if we 
replace A by A^ op in the results of the large-/3 approximation in Sec. 3. 
Similarly to the previous section, we can decompose Vjv(r) into four parts: 

V N (r) = V c (r) + B(N) +Cr + T>(r, N) + (terms that vanish as N -> 00), (23) 

where 

V c (r) = Im / dq a v (q), 24 

p r 7r Jci qr 

B(N) = Jim ^Re ^ dqe iqr {a v {q) - [a v (q)] N }, (25) 

C = ^ ( dqqa v (q), (26) 
2m Jc 2 

V(r,N) = V N (r) - [V c (r) + B(N) + C r\. (27) 

The integral contours C\ and C 2 are displayed in Figs. 2(i),(ii), respectively. 1 " 

The coefficient of the linear potential can be expressed analytically for (a),(b),(c). In the 
first two cases, the expressions read 

c <°> = ^ ( A S° P ) 2 . < 28 > 

cm = W (ASr) . (l + | r 7(1 + , „] _ (29) 



*This is the case when the number of quark flavors is less than 6 and all the quarks are massless. 

tWe conjecture that the expressions (24)-(27) are valid also beyond the 3-loop running, i.e. when the higher 
p n and a n are incorporated, as long as as = remains to be the IR fixed point when as(q) is evolved along 
Ci. 
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where 7(2;, r) = J T dtt x l e 1 represents the incomplete gamma function; Ai^ op and A^ op 
denote the Lambda parameters in the MS scheme; 5 = Pi/ Pi- In the case (c), C can be 
expressed in terms of confluent hypergeometric functions except for the coefficient of ci2, while 
the coefficient of can be expressed in terms of generalized confluent hypergeometic functions. 
Since, however, the expression is lengthy and not very illuminating, we do not present it here. 

The asymptotic behaviors of Vc(r) for r —>■ are same as those of Vqcd(^) i n the respective 
cases, as determined by RG equations. The asymptotic behaviors of Vc(r) for r —>■ 00 are given 
by — 47rCV/(/3or) [the first term of eq. (24)] in all the cases. 

As for B(N) and T>(r,N), we have not obtained simple expressions in the cases (b),(c), 
since analytic treatments are more difficult than in the case (a): we have not separated the 
divergent parts as N — » 00 nor obtained the asymptotic forms for r — > 0, r — > 00. Based on 
some analytic examinations, together with numerical examinations for N < 30, we conjecture 
that B(N) and T>(r, N) in the cases (b),(c) have behaviors similar to those in the case (a). 

Let us compare the "Coulomb-plus-linear" potential, Voir) +Cr, for the three cases when 
the number of quark flavors is zero. We also compare them with lattice calculations of the 
QCD potential in the quenched approximation. See Fig. 3. We take the input parameter 
for Voir) +Cr as a s (Q) = 0.2, which corresponds to A^°° p /Q = 0.057, A^°° p /Q = 0.13, 
A^ op /Q = 0.12.* Then, the scale for each lattice data set is fixed using the central value of the 

relation [18] A|^ op ro = 0.602(48), where tq is the Sommer scale. An arbitrary r-independent 
constant has been added to each potential and each lattice data set to facilitate the comparison. 
We see that Vc(r) + Cr for (a),(b),(c) agree well at small distances, whereas at large distances 
the potential becomes steeper as as(q) accelerates in the IR region, i.e. C^ a ' < < C^ c \ This 
feature is in accordance with the qualitative understanding within perturbative QCD [10,2,3]. 

*As well-known, when the strong coupling constant at some large scale, e.g. as(mb), is fixed, the values of 
^l-loop ^2-loop anc j ^3-loop jjjggj. substantially. As a result, if we take a common value of A-tjk as the input 

MS MS MS J MS " 

parameter, Vc(r) + Cr for (a),(b),(c) differ significantly at small distances, where the predictions are supposed 
to be more accurate. 
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The lattice data and Vc(r) + Cr also agree well at small distances, while they deviate at larger 
distances. More terms we include in Vc(r) + Cr, up to larger distances the potential agrees 
with the lattice data. If we increase the value of input as(Q), Vc(r) + Cr for (a),(b),(c) come 
closer to one another at r • A|^ op > 0.1. [The relation between Vc(r) + Cr for (c) and the 
lattice data remains unchanged.] 

5 Discussion 

In this section we discuss two issues: non- uniqueness of the decomposition of Vjv(r) and IR 
divergences of a n (n > 3). 

How to decompose V/v(r) for iV ^> 1 into the {r^ 1 , r° , r 1 , r 2 } terms is not unique. It is 
because Vjv(r) cannot be expanded in Laurent series about r = or r = oo due to logarithmic 
corrections. In fact, consider a function f(r) which behaves as const. + Kr for r and which 

Kr 2 

is Oil/r) for r ^> r±, where r\ represents a typical scale inherent in fir), e.g. fir) = — ; 

r + r\ 

then we may redefine Vc(r) = Vc(r) +f(r), C = C — K, V(r, N) = V(r, N) + Kr — f(r) + const. 
as the Coulomb part, the coefficient of the linear part, and the quadratic part, respectively. In 
particular, this redefinition changes the coefficient of the linear potential. 

On the other hand, we may consider the decomposition (23)-(27) to be an optimal decom- 
position for 10 < N < 100, on account of the following consideration. Suppose K is of the 
same order of magnitude as C. In the case r\ < A^, since Vc(r) + Cr is a good approxima- 
tion of Vjv(r) — B(N) for r < (see Fig. 1), Vc(r) + Cr cannot be a good approximation 
of Viv(r) — B(N) for ri < r < A^-. In the opposite case r\ > A^iL VMr) shows a linear- 

I VI o 1V1 o 

potential-like behavior for A^ < r ^ r\. Then it is not very appropriate to regard Vc(r) as 
the "Coulomb" part. 

As stated, a n for n > 3 contain IR divergences. In the computation of the heavy quarkonium 
spectrum based on potential-NRQCD formalism [19], IR divergences contained in Vqcd(^) are 
cancelled and the spectrum becomes finite at each order of the expansion in 0:5 [21,15]. Since IR 
divergences of Vqcd(^) originate from the separation of ultrasoft scale in the computation of the 
spectrum, it is natural to factorize the divergences from VQco( r ) by introducing a factorization 
scale fj,f (IR cutoff). In this case, VqcdM is rendered finite as well as dependent on /if. The 
IR divergences can be absorbed into a non-local gluon condensate. Thus, in the cases including 
a n for n > 3, it is sensible to investigate the truncated series Vjv(r) corresponding to Vqcd(^) 
regularized in this way. Another regularization scheme, which may be useful in comparison 
with lattice computations, is the resummation of a certain class of diagrams as done in [11], 
which turns the IR divergences into a finite contribution to the QCD potential. 

At the present stage, it is unclear to which regularization scheme or to which choice of \Xj 
in the factorization scheme the large-/3o approximation correspond. If we take the renormalon 
dominance picture rather strongly, we may expect that the regularization scheme dependence 
of the QCD potential is weak, and that the large-/? approximation makes sense quite generally. 
The full computation of a 3 will bring the status clearer on this point. 

One may consider that the leading ultrasoft logarithms of the QCD potential [22] are of the 
same order as the logarithms resummed by fa, so that they should be incorporated in the case 
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(c) of Sec. 4. It is achieved by replacing the corresponding V^-scheme coupling constant as 



in (24)-(26), where Ca = 3 is the second Casimir operator of the adjoint representation. We 
have checked that this contribution is very small and scarcely changes Vq{t) +Cr for the case 
(c) displayed in Fig. 3. This feature is consistent with the analysis of [5]. 

6 Conclusions 

We studied properties of the truncated perturbative series Vjv(r) of the QCD potential for 
N ^> 1; the perturbative expansion of Vqcd(^) is truncated at 0(ag) [N = 6n/(f3oas)}, where 
the term becomes minimal according to the estimate based on the NLO renormalon. V/v(r) was 
examined in the large-/3o approximation in Sec. 3. We decomposed Vn(t) into the {r _1 , r°, 
r 1 , r 2 } terms (with logarithmic corrections in the r^ 1 and r 2 terms) and analyzed properties of 
each term. The "Coulomb" and linear parts are finite as N — > oo, whereas the constant and 
quadratic parts diverge. We argued that the quadratic part can be interpreted as representing 
an C(A 3 r 2 ) uncertainty. For finite N, Vjv(r) is approximated well by the sum of the "Coulomb" 
and linear parts (up to a constant) for r < A" 1 and 10 < N < 100. In Sec. 4, higher-order 
terms of Vqcd(^) were estimated using the RG analysis. We decomposed Viv(r) into four parts 
and studied properties of the "Coulomb-plus-linear" potential. Analytic expressions for the 
linear potential are given in the 1-loop and 2-loop running cases. As we incorporate 1-, 2- 
and 3-loop running of ets(o), the linear potential becomes steeper, as well as the "Coulomb- 
plus-linear" potential agrees with the lattice data up to larger distances; cf Fig. 3. It is an 
interesting question whether the "Coulomb-plus-linear" potential converges toward the lattice 
data beyond 3-loop running. 

The linear potential is proportional to A|^r (as it should be unless it is zero, since there is 
no other dimensionful parameter). Naively one may think that such a linear potential cannot be 
produced within perturbation theory, since the expansion of A^g in «s vanishes to all orders. 
As eqs. (5), (6) show, however, the linear potential is indeed inherent even in a finite-order 
perturbative expansion of Vqcd(^), due to dimensional transmutation. In this regard, we note 
again that for certain finte N, Vjy(r) is approximated fairly well by the "Coulomb-plus-linear" 
potential (Fig. 1). 
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